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Abstract We use a composite gravitational galactic model consisting of a disk, a halo, 
a massive nucleus and a strong nuclear bar, in order to study the connections between 
global and local parameters in a realistic dynamical system. The local model is con- 
structed from a two-dimensional perturbed harmonic oscillator and can be derived by 
expanding the global model in the vicinity of the central stable Lagrange equilibrium 
point. The frequencies of oscillations are not arbitrary, but they are connected with 
all the parameters involved with the global model. Moreover, the value of the local 
energy is also connected with the value of the global energy. Low and high energy 
stars in the global model display chaotic motion. Comparison with previous research 
reveals that the presence of the massive nucleus is responsible for the chaotic motion 
of the low energy stars. In the local motion, the low energy stars show interesting res- 
onance phenomena, but the chaotic motion, if any, is negligible. On the contrary, the 
high energy stars do not show bounded motion in the local model. This is an indication 
of particular activity near the center of galaxies possessing massive nuclei. 
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1 INTRODUCTION 

In an earlier work (ICaranicolasll2002l - hereafter called Paper PI), we studied the connections be- 
tween the global and the local parameters in a barred galactic model. In the present research, we add 
a potential of a spherically symmetric nucleus and thus the total potential becomes 



^ r2 + « 2 Jr 2 [1 + (6 2 - 1) sin 2 <P] + c 2 



| In [r 2 [1 + (/3 2 - 1) sin 2 cf>] + c 2 ] - , (1) 



\Jr 2 + cl ' 



where (r, (f) are the usual polar coordinates. Equation ([T) describes the motion of stars in a barred 
galaxy with a disk, a halo, a massive nucleus and a strong nuclear bar (see below). Here Mj, Mb 
and M n are the masses of the disk, the bar and the nucleus respectively, while a,Cb,Ch and c n 
represent the scale lengths of the disk, the bar, the halo and the nucleus respectively. The strength 
of the nuclear bar is represented by the parameter b (b > 1), while the flattening parameter of the 
halo is represented by the parameter /3. Moreover, the parameter vq is used for the consistency of the 
galactic units. 
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In the system of galactic units used in this article, the unit of length is 1 kpc, the unit of time is 
0.97746 x 10 7 yr and the unit of mass is 2.325 x 10 7 Mq. The velocity and the angular velocity 
units are 10 km s _1 and 10 km s _1 kpc -1 respectively, while G is equal to unity. Our test particle 
is a star of unit mass (m = 1). Therefore, the energy unit (per unit mass) is 100 km 2 s -2 . In these 
units the values of the parameters involved are: a = 8,/3 = 1.3, 6 = 2, vq = 15, Aid = 9500, Mb = 
3000, M D = 400, Cb = 1.5, c n = 0.25 and — 8.5. The values of the above dynamical parameters 
remain constant during this research. 

We shall consider the case when the bar rotates clockwise, at a constant angular velocity Sib- The 
corresponding Hamiltonian, which is the well known Jacobi integral, in the rectangular Cartesian 
coordinates (x, y) is 

H 3 = \{pI+ Pi) + V) - \ul {x 2 + y 2 ) 

= Ej, (2) 
where p x and p y are the momenta per unit mass, conjugate to x and y respectively and 
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is the effective potential, while Ej is the numerical value of the Jacobi integral. If we expand the 
effective potential (0 in a Taylor series near the center, we shall obtain a potential describing local 
motion. 

The motivation of the present work is twofold: (i) to investigate the properties of global and 
local motion in the corresponding potentials. In particular, we shall express the coefficients of the 
local potential in terms of the global physical quantities entering the potential Q. A connection 
between the values of the global and the local energies will also be presented, (ii) to compare our 
numerical results with those obtained in Paper PI, where we only had the nuclear bar, while the 
massive nucleus was absent. 

The present paper is organized as follows: In Section 2 we study the properties of motion in the 
global model. The local potential, the connection between the local and the global parameters and 
the properties of the local motion, are presented in Section 3. We close with a discussion and the 
conclusions of this research, which are given in Section 4. 

2 PROPERTIES OF MOTION IN THE GLOBAL MODEL 

In Figure [TJ we can see the contours of the constant effective potential ([3J. The value of Sib is 1-25 
in the above mentioned galactic units. This value corresponds to 12.5 km s _1 kpc -1 . One observes 
that there are five stationary points, labeled L\ to L$, at which 

-5^ = , = . (4) 

ox ay 

These points are called Lagrange points. The central stationary point L\ is the minimum of $ c ff. 
At the other four points L2,L^, L4 and Lg, it is possible for the test particle to travel in a circular 
orbit while appearing to be stationary in the rotating frame. For this orbit, the centrifugal and the 
gravitational force precisely balance. The stationary points L4 and L5 on the x axis are saddle 
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Fig. 1 Contours of the constant effective potential (01. The values of all the parameters are 
given in the text. The five Lagrange points are labeled L\ to L5. 



points, while L2 and L3 are the maxima of the effective potential. The ann ulus bounded by the 
circle s through L2, L3 and L4, L5 is known as the "region of corotation" (see Binney & Tremainel 
l2008h . It is also important to note that the region of corotation in our dynamical system is located 
somewhere in the outer parts of the galaxy. 

We shall now proceed to study the properties of motion in the potential (|3j. We solve the equa- 
tions of motion in the (x, y) plane for the potential <[3j in which the nuclear bar rotates independently 
around the z axis. Assuming a clockwise rotation with a pattern angular velocity Ob, one can write 
the equations of motion in the form 



eft 



2(n 



b x r) + Ifibl r. 



(5) 



Decomposing into its x and y parts, we obtain 

dx 



u 



dy 



2n h y + nix , 

2n h ± + n 2 h y , 



(6) 



where the dot indicates derivative with respect to the time. 

All the numerical calculations are based on the numerical integration of the equations of mo- 
tion (O, which was made using a Bulirsh-Stoer routine in Fortran 95, with double precision in all 
subroutines. The accuracy of the calculations was checked by the constancy of the Jacobi integral (fJJ, 
which was conserved up to the eighteenth significant figure. 

In order to visualize the properties of motion, we shall use the classical method of the x — p x , 
y = 0, p y > Poincare phase plane of the Hamiltonian @. Figure|2]shows the structure of this phase 
plane for the global Hamiltonian ©, when Ej = —570. The particular value of the energy Ej was 
chosen so that in the phase plane x max — 10. There are regular orbits, forming the "right" retrograde 
set of invariant curves, as well as a triple set of islands. Furthermore, one can identify several sets of 
smaller secondary invariant curves embedded in the chaotic sea, which represent resonant orbits of 
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Fig. 2 The x — p x Poincare phase plane for the global Hamiltonian ©, when Ej = —570. 
The values of all the other parameters are as in Fig. Q] 

higher multiplicity. In addition to the regular region, there is a large unified chaotic domain formed 
by the chaotic orbits of the dynamical system. It is important to point out here that this chaotic area 
is obviously larger than that observed in the case in Paper PI , where the massive nucleus was absent. 
Also note that the presence of the massive nucleus gives rise to a variety of secondary resonances. If 
we set y = p y = in Equation ©, we obtain the limiting curve in the x — p x phase plane, which is 
the curve containing all the invariant curves for a given value of the Jacobi integral Ej. The limiting 
curve which is the outermost curve shown in Figure |2]is defined by the equation 

^pl + $ cff (x) = E 3 . (7) 

Figure |3a)-(f) shows six typical orbits of the global Hamiltonian system. The orbit shown in 
Figure [2 a) produces one of the "right" invariant curves. This orbit is a 1:1 resonant periodic orbit 
starting at the stable retrograde periodic point and it is nearly circular. Such orbits support the disk. 
The initial conditions are: xo = 3.89, yo = 0, p x o = 0. The orbit shown in Figure |3jb) produces 
the set of the three outer islands in the phase plane of Figure [2] This periodic orbit belongs to the 
family of the 1:3 resonant orbits. We observe that the shape of this orbit is elongated and therefore it 
supports the nuclear bar. The initial conditions are: xo = 9.32, yo — 0, p x o = 0. The orbit depicted 
in Figure|3jc) produces two of the three elongated islands of invariant curves which are embedded in 
the chaotic sea. This periodic orbit is characteristic of the 3:5 resonance and it has initial conditions: 
xq = —9.25, yo = 0, p x o = 0. In Figure|3ld) one can see a typical example of a periodic orbit which 
belongs to the family of 5:7 resonant orbits. The initial conditions are: xo = — 4.29, yo = 0,p x o = 0. 
The orbit shown in Figure|3le) is a complicated periodic orbit characteristic of the 5:9 resonance and 
produces a set of nine small islands of invariant curves inside the vast chaotic domain. This orbit has 
initial conditions: xo = —0.23, yo = 0, p x o = 73. It is evident that the last four types of periodic 
orbits and also the quasi-periodic orbits which belong to each family support both the disk and the 
bar structure of the galaxy. Finally, in Figure |3jf) we see a chaotic orbit with initial conditions: 
xo = —7.9, yo = 0, p x o = 0. The initial value of p y o was found in every case from the Jacobi 
integral (0. All orbits shown in Figure |3}a)-(f) were calculated for a time period of 100 time units. 
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Fig. 3 (a)-(f): Six typical orbits in the global Hamiltonian ©. The values of all the other 
parameters are as in Fig. |2] 



Let us now proceed in order to study the behavior of the orbits, in the global model, near the 
galactic center, that is close to the nuclear region. In order to visualize the nature of motion near 
the nuclear region of our dynamical system, we present Figure [4] which shows the x — p x , y = 0, 
p y > phase plane when Ej = —2700. Once more, the particular value of the energy Ej was 
chosen so that in the phase plane x max ~ 1, One observes a very interesting and complicated 
phase plane with regions of regular motion and a large, unified chaotic domain. In fact, there are 
two main regular regions consisting of invariant curves that are topological circles enclosing the 
two stable invariant points. The two stable invariant points D and R represent the direct (i.e. in the 
same direction as the rotation) and the retrograde periodic orbits respectively. These periodic orbits 
are characteristic of the 1:1 resonance and are similar to ellipses circulating around the origin. The 
rest of the regular region consists of several sets of smaller islands of invariant curves produced 
by quasi-periodic orbits belonging to resonances of higher multiplicity. The chaotic domain shown 
in the phase plane of Figure |4] consists of a large, unified chaotic sea and s ticky regions. It is well 
known that the phenomenon of stickiness is common in barred galaxies (see lCaranicolas & Karanisl 
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Fig. 4 The x — p x Poincare phase plane, when Ej = —2700. Details are given in 
the text. 

Emu aranis & Caranicolasll2002l) . The outermost curve shown in the phase plane of Figure |4] is 
the limiting curve defined by H(x,p x ) = Ej. 

Comparing the present results with those given in Paper PI, we understand that the presence 
of the massive nucleus dramatically changes the nature of motion in our barred galactic model. In 
Paper P 1 , the motion was regular and there was only one kind of invariant curve. The invariant curves 
were topological circles closing around the unique central invariant point. The corresponding orbits 
were box orbits. The outer invariant curves belonged to elongated box orbits that support the bar, 
but as we approached the central invariant point, the box orbits became more rectangular (see fig. 4 
of Paper PI for more details). On the other hand, in the present case, we have resonant periodic or 
quasi-periodic orbits producing several sets of islands of invariant curves. Furthermore, the majority 
of the phase plane shown in Figure|4]is covered by chaotic orbits. We will return to this point later 
in the discussion. 

Figure [5ja)— (j) shows ten typical orbits of the Hamiltonian system, when Ej = —2700. 
Figure |5ja) shows a periodic orbit starting at the position of the direct stable invariant point. This 
orbit produces an invariant curve which belongs to the main set of invariant curves and encloses the 
fixed direct invariant point. The initial conditions are: xq = —0.363, yo = 0, p x o = 0. In Figure|3Jb) 
we see a typical example of a 1:2 resonant periodic orbit, which produces the set of the invariant 
curves close to the limiting curve. This orbit has initial conditions: xq = 0.51, yo = 0, p x o = 32.8. 
The orbit shown in Figure |2c) produces the set of two small islands of invariant curves embedded 
near the main set of invariant curves around the retrograde stable point. This orbit represents the 2:2 
resonance and has initial conditions: xq = 0.152, yo = 0, p x o = 0. The orbit depicted in Figure|5Jd) 
produces one of the islands of invariant curves that intersects the p x axis. The initial conditions are: 
xq = 0.018, yo = 0, p x o = 38.6. In Figure|5je) we see a 3:4 resonant orbit with initial conditions: 
.To = 0.0756, yo = 0, p x o = 24.5. Figure|5jf) shows a periodic orbit characteristic of the 3:5 res- 
onance and produces a set of five small islands of invariant curves embedded in the chaotic region. 
This orbit has initial conditions: xq = 0.952, yo = 0, p x o = 0. In Figure |3g) one can observe a 
typical 5:5 resonant periodic orbit. This orbit produces the set of the five small islands around the 
invariant curves of the direct point and has initial conditions: xq = —0.609, yo = 0, p x o = 0. 
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Fig. 5 (a)-(j): Ten typical orbits in the global Hamiltonian (O when Ej = 
-2700. 



Order and Chaos in a Galactic Model with a Strong Nuclear Bar 



507 



Figure|3h) shows a periodic orbit which belongs to the 5:7 resonance family and produces a set of 
seven small islands of invariant curves. The initial conditions are: xq = 0.797, yo = 0,p X Q = 0. 
In Figure |5ji) we see a complicated periodic orbit characteristic of the 5:9 resonance, which pro- 
duces a set of nine tiny islands embedded in the chaotic sea. The initial conditions of this orbit are: 
xq = 1.025,yo = 0,p x o = 0. Finally, in Figure|5lj) we see a chaotic orbit with initial conditions: 
x n = 0.71,yo = QiPxO = 0. The initial value of p y0 was found in every case from the Jacobi 
integral (0. All orbits shown in Figure [5J a ) — CJ) were calculated for a time period of 100 time units. 

We have to point out that the nature of the local motion near the central region of the galaxy 
is much more complicated than the global motion. This can be justified by the fact that the phase 
plane shown in Figure |4] which corresponds to the local motion contains more resonant cases than 
the phase plane of the global motion, which is shown in Figure [2] Therefore, we may conclude that 
the presence of a massive nucleus near the galactic center gives rise to a large variety of resonant 
orbits of higher multiplicity. 

3 PROPERTIES OF MOTION IN THE LOCAL MODEL 



The local potential can be found by expanding the effective potential (01 in a Taylor series near the 
central stable Lagrange point L\, which coincides with the origin. Following this procedure and 
keeping only terms up to the fourth degree in the variables, we obtain the local effective potential 
which is 



U eB (Ax, Ay) = U cS (0, 0) + i [a (Axf + B (Ay) 2 



1 
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where we have set 



U, 
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in order to avoid large numbers. Setting for convenience: x 
U e s(0, 0), Equation <[8j becomes 
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As one can see from Equation (fTTT >. the coefficients of the local effective potential are functions of 
the physical quantities entering the global effective potential (01. 
The local Hamiltonian is 

H L = ±(pl+p 2 y )+V eS (x,y) = h h , (12) 

where p x and p y are the local momenta per unit mass, conjugate to x and y respectively, while hi, is 
the numerical value of the local energy. 

In order to connect the value of the global energy Ej with the value of the local energy hi,, 
we proceed as follows. The equation Ejq = ^ > eff(0,0) defines a point in the (x,y) plane, while 
Ej = $ c ff (x, y) defines a curve in the same plane. The global motion takes place inside this curve, 
which is known as the zero velocity curve. At the same time, /ilo = V c s (0, 0) defines a point in the 
(x, y) plane, while = U c h(x, y) defines a curve inside which the local motion takes place. This 
second curve is the local zero velocity curve. We only consider bounded motion, which means the 
zero velocity curves are always closed. The local energy is connected to the global energy Ej 
through the relation 

9 

ryC 

K = U eS (x, y) - U eS (0, 0) = — £ [$ cff (x, y) - $ cff (0, 0)] 

= ^(E 3 -E J0 ), (13) 



where 



a c b c n 



Let us now study the properties of the local motion. For the adopted values of the global pa- 
rameters and for fi b = 1.25, we find that: A = 1.4, B = 1.5, ai = 32.8, a 2 = 32.4, a 3 = 
32.8 and Oo = 0.009. For the value of the global energy Ej = —2700 we find, by using 
Equation Qj), that hi = 0.08452. It is amazing that this value of the local energy is much larger 
than the energy of escape for the corresponding local potential (see ICaranicolas & Varvoglisl[l984t 
ICaranicolas & Karanisll 1998b . which is given by 

(B — ill) 2 

h esc = y - 0) ■ (15) 

oai 

For the above values of the parameters, Equation dl~5t gives the value h csc = 0.01196. This value 
is less than hh/1. Therefore, one may conclude that near the massive nucleus we do not have local 
motion, or equivalently local orbits escape the nuclear region, because they possess a high value of 
local energy. In order to obtain an idea regarding the nature of the local motion near the nucleus, we 
must go to very low values of energies, that is when h <C hi,, where h is the value of the energy. 

Figure [6] shows the x — p x , y = 0, p y > Poincare phase plane when h = 0.012. The motion 
is regular and the phase plane has all the characteristics of the 1 : 1 resonance. There are two stable 
invariant points marked as D and R corresponding to direct and retrograde resonant periodic orbits 
respectively. The unstable periodic point gives an orbit, which in the absence of rotation is the p x 
axis. Thus, we observe that the local motion near the nucleus consists of a low energy 1:1 resonant 
periodic orbit. Moreover, as the entire phase plane is covered with invariant curves corresponding 
only to regular orbits, we conclude that near the nucleus we do not observe any kind of local chaotic 
motion. This situation is completely different from that displayed in Paper PI, where the nucleus 
was not present. In that case, all orbits in the local potential were box orbits, without any resonance 
phenomena. 
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Fig. 6 The x — p x Poincare phase plane for the local Hamiltonian (fTSI i. when h = 0.012. 




x x 

(c) (d) 

Fig. 7 (a)-(d): Four typical orbits in the local Hamiltonian (fT2l , when 
h = 0.012. 

Figure [7ja)— (d) shows four orbits of the local Hamiltonian (TlZb . Figure [Vj a) and (b) show two 
periodic orbits starting at the direct and the retrograde periodic points respectively. The initial con- 
ditions for the orbit shown in Figure [TJ a) are: xq = —0.0845, yo — 0, p x o = 0, while for the orbit 
shown in Figure|3b) are: xq = 0.0885, yo = 0, p X Q = 0. In FigureQc) we observe an orbit starting 
at the unstable periodic point with initial conditions: xq = —0.0119, yo = 0, p x o = 0. Finally, 
in Figure ITJd) we see a box orbit which produces one of the outer invariant curves shown in the 
phase plane of Figure [6] The initial conditions for this orbit are: xq = 0.1425, yo = 0, p x o = 0. 
The initial value of p y o was found in every case from the energy integral (TT2l . All orbits shown in 
Figure |7ta)-(d) were calculated for a time period of 100 time units. 
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Fig. 8 Time evolution of the L.C.E. for two chaotic orbits. The curve labeled G corre- 
sponds to an orbit in the chaotic domain of the phase plane of Fig. [2] while the curve 
labeled L corresponds to an orbit in the chaotic domain of the phase plane of Fig. |4] 



We must clarify to the reader that the Taylor expansion (TTOt is only valid when 



« 1, 



6V 



«1, 



x 2 + y 2 



«1, 



< 1. 



(16) 



It is also important to note that, for a given value of the global energy Ej, a corresponding value of 
the local energy can be obtained through relation dT~3T >. It is obvious that this value of the local 
energy does not have any physical meaning if all the relations ( TToT l are not satisfied. 

In order to better estimate the degree of chaos displayed by the chaotic orbits shown in Figures [2] 
and [4] we decided to compute the m aximum Lyapunov Characteristic Exponent (L.C.E.) (for details 
see lLichtenberg & Liebermanl 19921) . The results are shown in Figure[8] The curve labeled G shows 
the evolution of the L.C.E. for a chaotic orbit with initial conditions in the chaotic region of Figure|2] 
The particular values of the initial conditions of this orbit are as in Figure |3jf). On the other hand, 
the curve labeled L shows the evolution of the L.C.E. for a chaotic orbit with initial conditions in 
the chaotic region of Figure |4] The particular values of the initial conditions of this orbit are as in 
Figure|51j). We observe that the value of the L.C.E. corresponding to the local system of Figure|4]is 
about ten times the value of the L.C.E. corresponding to the global system of Figure [2] Th erefore, 
one can say that in this case we have not only fast chaos (see ICaranicolas & Vozikisll 19871) . where 
the L.C.E. was on the order of unity, but also very fast chaos where the L.C.E. is about three times 
larger. It is evident that this is an indication of strong nuclear activity near the vicinity of the galactic 
center. 



4 DISCUSSION AND CONCLUSIONS 



One of the most important approaches in order to understand and reveal the dynamical behavior 
of a galactic system is based on the knowledge of the chaotic versus ordered nature of orbits. 
In the present research, we tried to achieve this using a global and a local potential describing a 
barred galaxy with a massive nucleus, a disk and prolate halo components. The corresponding lo- 
cal potential was found by expanding the global potential around the central stable Lagrange point, 
which coincides with the origin, in a Taylor series and keeping only terms up to the fourth degree 
in the variables. This local potential is a potential made up of a two-dimensional perturbed har- 
monic oscillator. The study of the motion in those potentials has been an active field of research 
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over t he last decades (see for examples ISaito & IcMmuralll979t ll nnane r] ll985UCaranicolaslll984 
11994 l2000t ICaranicolas & Karanislll999l ICaranicolas & Vozikisll2002l) . In later years, the study of 
the properties of motio n in those dynamical systems has been e xplored in detail, using precise and 
modern analytical (see lElipell2000l 1200 it lElipe & Depritll 19991) or numerical (see lLara et al.lH999t 
iKaranis & Caran icolas 2002) methods. 

As expected, the local parameters and the corresponding local energy are functions of all the 
involved global parameters and the global energy. Our numerical experiments, in the global dynami- 
cal system for high values of the energy, suggest that more than 62% of the tested orbits are chaotic. 
Comparing the outcomes with those of Paper PI we see that, in the present case, we have a sharp 
increase of the total chaotic motion. It is evident that the chaotic domain in Figure [2] is much more 
extended. This means that the additional nucleus affects not only the region near the center of the 
galaxy, but also areas far from it. Another interesting observation is that the presence of the nucleus 
does not seem to affect the figure-eight quasi-periodic orbits, which are the building blocks of the 
barred structure of the galaxy. This result agrees with recent observations made by the Hubble Space 
Tele scope, which revealed that a number of Seyfert galaxies display strong nuclear barred structure 
(see lRegan & Mulchaevl[l999l) . On the other hand, the global motion near the massive nucleus has 
the characteristic of a 1:1 resonance. Furthermore, the chaotic domain seems to be larger than that 
observed in Figure [2] This is natural because of the presence of the nearby massive nucleus. All the 
above results are completely different from the corresponding results obtained in Paper PI, where 
the motion was regular and all the orbits in the global system near the galactic center were box orbits. 

The effect of a central mass concentration (CMC) on the phase space of global galactic models 
has already been st udied extensively i n seve ral earlier research works (see Caranicolas & Zotos 
2011; Zotos 2012). lHasan & Normanl dl99dt) have shown that a black hole or a CMC in a barred 
galaxy can dissolve the barred structure. Moreover, the fundamental orbits of the system (called 
B orbits), supporting the bar, are orbits elongated in the direction of the bar. As the black hole's 
mass increases, an inner Lindblad resonance (ILR) appears and moves outward. Thus, the B orbits 
disappear when the ILR reaches the end of the bar. This effect is significantly increased if the bar 
is thinner. Orbits that stay closer to the central mass, with smaller Jacobi co nstants, can change 
their character from regular to chaotic more easily. Furthermore, the study of lHasan et al l (1 19931) 
indicates that, as the mass of the CMC increases, the stable regular orbits in the region where the bar 
potential competes in strength with the central mass potential, especially aroun d the region of the 
ILR, w ill no longer be present and their position will be occupied by chaotic orbits. IShen & Sellwoodl 
(12004 have conducted a systematic study of the effects of CMC on bars, using high quality TV-body 
simulations. They have experimented with both strong and weak initial bars and a wide range of 
the physical parameters of the CMC, such as the final mass, the scale length and the mass growth 
time. Their outcomes suggest that, for a given mass, compact CMCs (such as super-massive black 
holes) are more destructive to barred structures than are more diffuse ones (such as molecular gas 
clouds in many galactic centers). They have shown that the former are more efficient scatterers 
of bars supporting regular elongated orbits, that pass close to the galactic center and, therefore, 
decrease the percentage of the regular orbits and increase the area of the chaotic region in the phase 
space. All the above research outcomes are in agreement with the present findings, which strongly 
indicate that the massive nucleus is a very important parameter of the dynamical system, since it is 
responsible not only for the chaotic motion of low energy stars but also for the several resonant orbits 
of higher multiplicity that appear in the local model. In particular, Figures[2]and[4]in the current pape r 
remarkably resemble the Poincare surface of section plots in figure 1 1 of Shen & Sellwoodl d2004 . 
where a realistic self-consistent barred model was adopted. 

Let us now consider the local motion. The local potential (TlOb has all the c haracteristics of the 
1:1 resonance. Such potentials are known as perturbed elliptic oscillators (see |Pepritlll99ll) . It is 
very interesting that for the corresponding local energy given by Equation ( [13) , the local motion 
is not bounded. In order to obtain bounded motion, one must approach very low values of local 
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energies. Strictly speaking, the presence of the massive nucleus makes the local energy increase 
dramatically, so that considering all the zero velocity curves results in unbounded local motion. For 
values of local energies slightly lower than the value of the energy of escape, we observe resonance 
phenomena, namely the 1:1 resonant periodic orbits. Extensive numerical calculations in the local 
potential suggest that chaotic motion was not observed. 

Bearing all the above in mind, we can say that this situation is entirely different from that inves- 
tigated in Paper P 1 , where the properties of local motion were the same as those of the global motion 
near the galactic center (see figs. 4 and 5 in Paper PI). In our case, the corresponding local motion 
does not appear to exist and only very low energy local motion seems to be present. Furthermore, 
the fact that the corresponding value of the L.C.E., which was computed in the chaotic domain near 
the nuclear region, is much larger than unity suggests that there is a particular and strong activity in 
the central parts of galaxies with strong nuclear bars. 

We consider this paper to be an initial effort in order to explore and reveal the dynamical struc- 
ture of the system in more detail. As the present results are positive, further investigation will be 
initiated to study all the available phase space. Moreover, our gravitational model will be suitably 
modified, in order to be able to describe the properties of motion in a Hamiltonian galactic system 
of three degrees of freedom. 
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